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ABSTRACT 

There is growing evidence that the estimations of the beaming Doppler factor 
in TeV BL Lac object based on the Self Synchrotron Compton (SSC) models are 
in strong disagreement with those deduced from the unification models between 
blazars and radio galaxies. When corrected from extragalactic absorption by the 
diffuse infrared background (DIrB), the SSC one-zone models require very high 
Lorentz factor (around 50) to avoid strong 7 — 7 absorption. However, the statis- 
tics on beamed vs. unbeamed objects, as well as the luminosity contrast, favors 
much lower Lorentz factor of the order of 3. In this paper, we show that for the 
special case of Markarian 501, the need for very high Lorentz factor is unavoid- 
able for all one- zone models where all photons are assumed to be produced at the 
same location at the same time. Models assuming a double structure with two 
different beaming patterns can partially solve the problem of luminosity contrast, 
but we point out that they are inconsistent with the statistics on the number of 
detected TeV sources. The only way to solve the issue is to consider inhomoge- 
neous models, where low energy and high energy photons are not produced at 
the same place, allowing for much smaller Lorentz factors. It implies that the jet 
is stratified, but also that the particle energy distribution is close to a monoen- 
ergetic one, and that pair production is likely to be significant. The implications 
on relativistic jet physics and particle acceleration mechanism are discussed. 
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Subject headings: radiation mechanisms: non-thermal — gamma rays: theory — 
galaxies: active — galaxies: BL Lacertae objects: general — galaxies: jets 

1. Introduction 

It is now admitted that the blazar phenomenon is due to relativistic Doppler beaming 
of the non-thermal jet emission taking place in radio-loud Active Galactic Nuclei (AGN) 
where the jet axis is closely ahgned with the observer's line of sight. They exhibit an 
important level of optical polarization, a fiat radio spectrum, a strong variability in all 
frequency bands and a very broad spectral energy distribution (SED) ranging from the radio 
to the extreme gamma ray band. The SED consists typically in two broad components. 
In the so-called Synchrotron-Self-Compton process (SSC) model, the lowest energy hump is 
attributed to the synchrotron emission from relativistic electrons and/or positrons, and the 
highest one is attributed to the Inverse Compton (IC) process of the same charged particles 
onto the synchrotron photons and/or external photons. The blazar class of objects includes 
both the Flat Spectrum Radio Quasars (FSRQ) and the BL Lac sources (or Lacertids), 
depending respectively on the existence or the lack of detectable emission lines in their 
spectra. Following Chiaberge et al. (2000), one can define two classes of BL Lac objects 
(which are most probably two extreme cases in a continuous distribution) : the LBL or 
"red" BL Lac, for which the synchrotron component peaks in far IR to optical, and the IC 
component peaks in the MeV-GeV range, and the HBL or "blue" BL Lac, for which the 
synchrotron component peaks in the UV-X range, and the IC component peaks above 10 
GeV. The most extreme objects up to now are those whose non thermal emission extends 
up to the TeV range, the so-called TeV blazars. The two main prototypes are Mrk 421 
(Punch et al. 1992) and Mrk 501 (Quinn et al. 1996), two radio-loud AGN relatively close 
to us and roughly at the same distance, Zg ~ 0.031 and Zg ^ 0.034 respectively. Five other 
TeV detections have been repeatedly detected (lES 1959-1-650, PKS 2155-304 (Aharonian et 
al. 2005a), lES 1426+428 and PKS 2005-489 (Aharonian et al. 2005b), and lES 2344+514 
(Aharonian et al. 2004)). All of them are Lacertids, although it is not clear up to now 
whether only BL Lac objects do emit TeV radiation or if this is due to a selection effect. As 
a matter of fact, BL Lac objects appear to be much more numerous than quasars and the 
closest blazars all belong to this class. A high sensitivity threshold would strongly bias the 
detection toward the closest sources. Furthermore it is well-known that TeV photons are 
absorbed by the Diffuse Infra Red Background (DIrB) to create electron-positron pairs, and 
it is not obvious whether even the closest quasar, 3C 273 {zg ~ 0.158) could be detected in 
the TeV range. 
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One-zone SSC models assume that highly relativistic particles are injected in a spherical 
zone, where they cool by emitting synchrotron radiation and by Inverse Compton process. 
The models require specifying the source radius, the magnetic field, as well as the density 
and the energy distribution of the emitting relativistic particles. The latter is most often 
assumed to be a power law or a broken power law {e.g. Marscher 1983; Tavecchio et al. 
1998). It turns out however that the computation of emitted radiation is not compatible 
with the hypothesis of a static source, because in most cases the photon density would be so 
high that all TeV photons should be absorbed to form electron-positron pairs. Furthermore 
the time variability is so short (down to 15 minutes in some cases, Gaidos et al. (1996)) that 
it is incompatible with a spherical static source through the causality argument. This leads 
to assume that the source is moving with a relativistic bulk velocity v = (3c. The effect of 
relativistic bulk motion is entirely described by the Doppler beaming factor 5 = l/r(l — /3/x), 
where F = (1 — /3^)~-^/^ is the usual Lorentz factor and fi = cos 9 is the cosine angle of the jet 
according to the observer's line of sight. The Doppler effect shifts all frequencies by a factor 
6 and all specific intensities by a factor 6'^. So the actual photon density in the jet frame is 
much lower than what would be deduced for a static source. The relativistic motion has been 
invoked for a long time (Rees 1966) to solve a similar issue for radio emission of quasars. 
Namely the brightness temperature is so high that, for a static sources, the relativistic lep- 
tons emitting synchrotron radiation should have cooled immediately through the so-called 
"Inverse Compton catastrophe" (Rees & Simon 1968). Again the relativistic motion can fix 
this issue, because the actual photon density in the jet frame is much lower when taking into 
account the Doppler amplification. This beautiful theoretical explanation has been later 
confirmed by the discovery of superluminal motion, which requires Lorentz factors at least 
as great as the observed apparent reduced velocity (3app = ^app/c (for a review, see Zensus 
1997). 

For > /3, corresponding to ^ < 1/F, one has 1^5^ 2F, whereas 5 ~ 1/F outside this 
interval. It means that for a few beamed Doppler sample of boosted sources, one expects a 
lot of unbcamcd and not amplified counterparts. It is natural to think that the unbcamcd 
counterparts of bright quasars are the weaker radio galaxies, whose jet is thought to make 
a larger angle with the line of sight. Particularly it has been proposed that the unbeamed 
counterparts of BL Lac object could be a sub class of radio galaxies, the so-called Fanaroff- 
Riley I (FRI) radio galaxies (Urry & Padovani 1995). These are characterized by a rather 
faint, weakly beamed, and core-brightened radio jet. Statistical studies of radio and X-ray 
AGN samples have confirmed the possibility of such an association. The inferred beaming 
factors seem to imply a rather modest value of the bulk Lorentz factor, of about 3. However, 
the modeling of SSC radiation by one zone models requires much higher values : following 
the authors, they range from 10 to 50 (Tavecchio et al. 1998; Konopelko et al. 2003; Sauge 
& Henri 2004). The highest value seems to be needed when one takes properly into account 
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the extragalactic absorption. The problem is further comphcated by the absence of clear 
superluminal motion in TeV blazars, together with a rather modest brightness temperature, 
which implies also a low Lorentz/Doppler factor (Edwards & Finer 2002; Finer & Edwards 
2004). All these contradictory facts lead to what we call here the "Bulk Lorentz factor crisis 
of TeV blazars" . 

The aim of this paper is first to ascertain this crisis. Wc will first show that all one- zone 
SSC models imply high Lorentz factors, only with the argument of 7 — 7 absorption and 
discarding any variability argument. Then wc will recall the arguments for low Lorentz fac- 
tor, based on general geometric properties of the Doppler boosting. We will show that the 
explanations based on two different structures, with a possible deceleration of a fast spine 
responsible for TeV emission, are not satisfactory concerning the statistics of TeV blazars. 
We argue that the best solution is to admit the low Lorentz factor constraint, abandoning 
the one-zone assumption. We will show that this conclusion has important consequences 
regarding the jet physics and the particle acceleration mechanism. 



2. The case for high Lorentz factor 

In the following, we will develop the need for high Lorentz factors for one-zone models, 
with the fewest theoretical assumptions and relying only on observational data. We will only 
assume that the SSC process is at work, with the usual assumptions of one-zone models : the 
rclativistic particles are assumed to be injected in a spherical homogeneous "blob" of radius 
R, moving at a rclativistic velocity characterized by the Lorentz factor F and a corresponding 
Doppler factor 6. The blob is filled with a tangled magnetic field of constant strengh B. We 
will refer to all quantities expressed in source rest frame by a star and quantities in observer's 
frame are not labeled. All energies are expressed in reduced unit of mgC^. Throughout this 
paper, we express the Hubble parameter by Hq — 100/ikms~^ Mpc~^ and assuming h to be 
equal to h — 0.65. 



2.1. The synchrotron and IC differential Luminosity 

Inspection of the TeV blazars spectra shows that the IC spectra reaches their maximum 
luminosity at some peak energy Sc, which is of the order of 10^ for TeV photons. This energy 
corresponds to an energy e* = ecS~^ in the blob frame. We will consider only the particles 
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emitting this typical energy via the IC process, which have a typical individual Lorentz factor 
(in the blob frame) 7c, which must be greater than e*. 

We then define another typical energy Eg, that is emitted by synchrotron process by the same 
particles. It can be expressed in the blob frame as e* = (B/Bq)'-^'^ where Bq = 3Bc/2 and 
Be = 2TTmlc^/eh ^ 4.41 x lO^^G is the usual "QED critical magnetic field strength". One 
has then Eg = 6{B/ Bo)j'^. Synchrotron spectra of TeV blazars are typically peaking in the 
1-100 keV range so that ~ 10"^ - lO-\ 

Synchrotron photons are up-scattered at high-energy via Inverse Compton process. It has 
been stressed by various authors that, giving the observed energies of IC and synchrotron 
photons, the collisions between the most energetic particles and the peak synchrotron photons 
take place in the Klein- Nishina regime, that is £*7c ^1- In this condition, the particle 
(electron or positron) gives all of his energy in a single interaction. It follows that 7c ~ £*. 
This gives an estimate of the magnetic field strength, 

B = Boj^, = 6B,% (1) 

which is only valid in the Klein-Nishina scattering regime. 

The differential synchrotron luminosity L^ g per unit reduced energy emitted by a pop- 
ulation of particles of energy 7 at the energy Sg with differential energy number of particles 
dA^/d7 reads 

L^sg) = —isg) = 6 ^{Sg) = S (2) 

The total power lost per particle of energy 7 is given by the well-known relation dL*/dA^ = 
(4/3)c(TTh7^W/s, and we obtain, 

4 £^ dN 

L„fe) = 53«.TKH'B-^ (3) 

where Wb = B^/8tt is the usual magnetic energy density. Combining with equation (1), we 
can write : 

LsA^s) = S'^ca^^Bf-^^. (4) 

OTT Ec d7 

We can compute the differential IC luminosity L^jci^c) in the same way using an expression- 
similar to Eq. 2, but replacing the magnetic energy density by the photon energy density. 
However wc have to take into account that the Klcin-Nishina cut-off reduces the effective 
energy density available for Inverse Compton scattering. We thus define a new characteristic 
energy, corresponding to the synchrotron photon energy at the limit between the Thomson 
and Klein-Nishina regime for particles with an energy 7c. This energy is — l/7c, i.e. in 
the observer's frame 

St = Sye,. (5) 
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Photons of this energy will also be the main contributors for absorbing photons of energy Sc 
to create electron/positron pairs. If we neglect the Klein-Nishina contribution above St, the 
total power lost per particle of energy 7 writes dL^/dN — (4/3)c(7Th7^M^ph*) where 



1 r A (R^ 



For a power- law spectrum uF^, oc u'^ with > 0, a simple calculation gives 

The coefficient f3 can be replaced by another numerical coefficient close to 1 as long as the 
vF^ spectrum is growing with energy. The differential IC luminosity reads then : 

LeMe.)^S-^p-'^s.L,M)^. (7) 

Comparing equation (4) and (7) we can now estimate the radius of source as a function 
of observed luminosities and the unknown Doppler factor : 

We will now use this radius estimate to compute the 7 — 7 optical depth for the photons of 
energy Ec- 



2.2. The 7 — 7 photon opacity 

As we mention above, gamma-ray photons of energy are mainly absorbed by photons 
of energy £*~^ = creating pairs. So the same soft photons control both the amount of IC 
process and the absorption of IC photons. The absorption probability (or opacity) per unit 
path length of a photon of energy £* ^ 1 due to pair production in the case of a power-law 
SED is given approximately by 

d 

where £j^^{e*) is the free mean path of the photon and n{e*) the differential photon density 
per unit of reduced photon energy e*. In the framework of one- zone model, the typical 
interaction scale is of the order of the size of source R. It follows that the typical 7 — 7 
optical depth writes, 

R 

^77(£*) ~ " a^jRaThsM^D- 

ijj[e ) 
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The function a^^ (Svensson 1987; Coppi & Blandford 1990) depends on the index (5 of the 
power-law of the spectral soft photon density expressed in vFy oc z/^ form. A commonly 
used value of a^^{l3) is 0.2 or 0.25. More precisely we have (Svensson 1987) 



a (3) = A^-^6x '^'^'-^K''-^^''-^^''-^^^ (9) 



The differential energy density number of particle is given as a function of the differential 
luminosity by 

so we get finally the optical depth as a soft compactness at the energy et- 

T^^{e ) a^^^p) 4T,rn,cm ^^^^^ Anm^c^R ' ^ ^ 

Using our estimate on the source radius R equation (8), we obtain 



r,,(s.) = r^a,,(/3)^^ig^ (12) 



1/2 

1/2 



/ L£^s{^t)Le,Ic{^ 



where we introduce the modified function a^^{l3) as Q;^-y(/3) = a^^\f^. Values of a^^ and 
for some /3 are tabulated in table 2. 



2.3. Constraints on the local synchrotron spectral shape 

Equation (2) shows that if we are able to measure the position in frequency and fiux 
of both the synchrotron and the IC peak, then we can evaluate the optical depth to 7 — 7 
absorption at the IC peak as a function of the assumed Doppler factor value. This optical 
depth is controlled by the synchrotron luminosity at the frequency St — jsc- We can 
use this relation either by assuming some Doppler factor and evaluate the optical depth, 
or constrain the value of 5 by limiting the value of r-y-y. We can define Tmax, the Compton 
dominance parameter, as the ratio of IC luminosity's peak to the synchrotron one, 

^ _ K^c(^c)]niax _ ej.,^ic{e^ 

'max rr-i/M r / \' V 

We can rewrite equation (12) to express the luminosity at = as a function of the 

optical depth and the rmax parameter. We finally obtain, 

r^^(£--) 6eh el 



-max 



Q;^-y(/3) (7Th?77.eC^/^ g-max 



(14) 
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Equations (5) and (14) can be considered as a system of two parametric equations of the 
curve giving s^Lg ^(^t) as a function of St- Eliminating 5 between the two previsously cited 
equation, one gets the following expression : 



(kh 



a/2, 



(15) 



For nearby sources, the luminosity distance writes di{z) ~ cz/Hq and previous expression can 
expressed in term on flux F instead luminosity L using the well known relation F — L/Aird'^, 



dFs 



= 3.4 X 10 



-36 r4 



h T, 



77 



'^max\ ^max 



z a. 



77 



(16) 



For an observed SED and a given value of the opacity parameter r^^, the only remaining 
unknown quantity in the previous equation is the beaming Doppler factor 6. Each value of 6 
gives a point in the log^^Q z/ — logj^Q(z/Fi,) plane lying on straight line of slope 2, the level of the 
curve depending only on the value of r^-y. Intersection of the synchrotron spectrum with the 
straight line directly constrains the minimum value SminiT-y-y) of the beaming Doppler factor 
required to avoid the 7 — 7 absorption with an opacity value of Tjj of the Inverse Compton 
bump (at the peak frequency). 



2.4. Application to Mctrkarian 501 

We apply this calculation to the case of the Mrk 501 object during the period of the 
1997 April 16 where the Beppo-SAX satellite (Plan et al. 1998) and the CAT imaging At- 
mospheric Cerenkov Telescope (Djannati-Atai et al. 1999; Barrau et al. 1998) have recorded 
simultaneous data (see figure 2). All observational parameters we need in the equation (16) 
are reported Table 1 . We consider the two cases where we take into account, or not, the 
attenuation of the high energy component by cosmic diffuse infrared background (DlrB). 
This effect consists in the interaction of emitted gamma rays during their travel through 
the Universe with the photon field of the diffuse infrared background (DlrB) to create pairs 
(Gould & Schreder 1967a,b; Stecker, de Jager, & Salamon 1992; Vassiliev 2000). The tail of 
the high energy spectra is then de-reddened using the method described in Sauge & Henri 
(2004). This situation changes the position of Inverse Compton peak and the Compton 
dominance parameter Tmax- In this case, for r^^ = 1, we obtain both in the reddened and 
the de-reddened case Smini^) ~ 50 (see figure 3). 

Given equation 16, the position of the line constraining 5 depends on the value of 
{e^^y/L^jc- It turns out that, also the IR un-folding of the spectrum changes both quan- 
tities, the previous ratio depends only slightly on the level of assumed absorption. The 
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value of 5inin(l) are thus quite similar in the two cases because when we correct the Inverse 
Compton bump, position of the maximum moves both in luminosity and in frequency. This 
effect could be clearly seen on figure 3, where the difference between the two panels is hardly 
perceptible. 

Note that in fact the level of the curve depends implicitly also on the value of modified 
power-law index of the spectrum (5 (see eq. [16]). In our case, we choose a value (5 — 0.5 
directly measured on the SED. 



Another constraint can be derived from the observation of short variability timescale. 
The classical argument is that a spherical static source cannot be variable on a timescale 
smaller than R/5c. So one gets an upper bound of radius of the source R ^ Rvar.imu = ^ctvar- 
Combining previous inequality to the equation (8) and expressing all the quantities in their 
fiducial units, we finally get a constraint similar to the one obtained in the previous section 
for the local synchrotron shape (see eq. [16]) : 



Taking a characteristic variability timescale of roughly 15 min, we obtain the left solid thick 
line displayed on the figure 3. It appears that this constraint is less restrictive than the 
previous. In context of homogeneous modeling, it gives a minimum value for Doppler factor 
of 6-8 and 8-10 for the reddened and the de-reddened case respectively. 



In this section we shortly review all the arguments and pieces of evidence in favour of 
moderate or low values of the bulk Lorentz factor. 



Observations at the VLBI scale (^mas) show that blazars often display superluminal 
apparent velocities. This phenomenon predicted by Rees (1966) is expected for relativistic 
moving sources which is highly beamed and closely aligned with the observer's line of sight. 



2.5. Constraint from the variability timescale 




(17) 



3. 



The case for low Lorentz factor 



3.1. Absence of superluminal motion at peirsec scale 
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For a component moving along the jet axis at a reduced speed (5 — v/c and making an angle 
6 from the line of sight, the apparent transverse velocity measured by the observer is : 

If > /^crit = a/2/2 and 9 is such that sin2^^ > (F^ — the motion will appear to be 

superluminal i.e. /Sapp > 1- Expressed in degrees the latter condition writes > 6*^^ = 
0.28 (F/10)-2 deg. 

As a matter of fact, VLBl/VLBA campaigns have not clearly succeeded in finding su- 
per himinal motion at the parsec scale for any TeV blazars (Edwards & Finer 2002; Finer & 
Edwards 2004). Observed radio components seem to be stationary or subluminal, requiring 
low or moderate values of the Lorentz factor (F 2 — 4). 

The absence of superluminal motion could be explained by a very close alignment of the jet 
with the hue of sight. Indeed, following the previous expressions, if ^ ^ 1/2F^ the apparent 
velocity is always smaller than c, and object appears to be subluminal despite the large 
value of F. But it this case, a simple statistical argument based on the density number of 
unbeamed counterparts rule out this possibility as we will see in the next section. 



Moreover, derived value of the brightness temperature of the VLBI core is in the order 
of 10^*^"^^ K and lie well below the usual Inverse Compton limit of ~ 10^^~^^ K necessary 
to avoid the " Inverse Compton catastrophe" , i. e. situation where ultra- relativistic particles 
suffer from dramatically Compton cooling in a very short time. Finer & Edwards (2004) have 
concluded that the jet should be only mildly relativistic at parsec scale. They propose that 
the TeV emitting inner jet is strongly decelerated before reaching the parsec scale. However 
we will see in the following that the existence of the highly relativistic motion is challenged 
by other observational facts concerning the statistics of beamed vs. unbeamed sources. 



3.2. Number of beamed sources 
in the BL Lac/FR-I unification pctradigm 

As we said in the introduction, the blazar phenomenon arises from a close alignment of 
jet axis with the observer's line of sight. Following this scheme, one expects the existence of 
sources sharing the same physical properties {i.e. intrinsically the same objects), but viewed 
at larger angle. It has been proposed that Fanaroff- Riley radio galaxies can be the unbeamed 
parent population of blazars and particularly, FR-I galaxies can be the counterparts of 
Lacertids (Urry & Fadovani 1995). 
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The unification hypothesis can be tested on samples of objects both by their luminosity 
ratio and by their spatial density. Doppler beaming effect enhances the intrinsic bolometric 
luminosity by a factor 5^. The Doppler factor itself varies from ^max = 2r for a jet pointing 
exactly toward the observer, to ^min = l/F for jets lying close to the sky plane. Although 
the exact definition of what is a beamed object can be somewhat subtle, one can estimate an 
order of magnitude of the number of such objects. It is easy to see that the solid angle for 
which the Doppler factor is larger than some given value 5q (where of course < 5q <2V 
is 

a^2.(i-,„).i^(i-±) (19) 

where we used ~ 1 + l/(2r^), and hence the fraction of sources with a Doppler boosting 
larger than 5q is approximately 

if we assume always two symmetrical jets. 

So any "beaming criterion" imposing a Doppler factor larger than some sizable fraction of 
r will give a fraction of beamed sources of the order of F"^. For one beamed source, one 
expects thus around F^ unbeamed sources. It turns out that careful statistical studies do 
indeed confirm the association of BL Lac objects with FR-I galaxies, but they converge 
toward a much lower Lorentz factor than what is expected from the gamma-ray emission. 
For X-ray selected BL Lacs, which comprise all known TeV blazars, the inferred density ratio 
is 1:7, corresponding to a bulk Lorentz factor around 3.5 (Urry & Padovani 1995). 

In the same way, we can examine the hypothesis that the lack of detection of super- 
luminal motion would be due to a close alignment of the jet with the line of sight. As we 

note in the previous section, a beamed source with 9 ^ 1/2F^ can not appear superluminal. 
The cone substained by this angle corresponds to the solid angle ^ 7r/4F'^. The ratio /' 
between the density of unbeamed sources and subluminal beamed sources one is thus given 
by 

/' = « 16Fl (21) 

'^TcV.sub 

Observations show that 5 TeV blazars do not clearly display superluminal motions for ^ 
0.047 (Edwards & Finer 2002; Finer & Edwards 2004). Assuming 5 = 2F = 50, it corresponds 
roughly to a volume of 0.043 Gpc^ and then to a density of subluminal TeV blazars of 
?^TeV,sub = 117 Gpc~^. Then, the density of expected unbeamed counterparts would be 
fiFK-i ~ 7.3 X 10^ Gpc~^ which is absolutely unreasonable. 



3.3. Luminosity ratio of beamed sources 
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in the BL Lac/FR-I unification pciradigm 



Another constraint can be derived from the luminosity ratio between the Lacertids 
sources and the FR-I ones. Supposing the same assumption as above (see previous section 
3.2), ie Lacertids and FR-I are on average the same intrinsic objects but viewed at different 
angles, the bolometric luminosity contrast between the two parents population (Lacertids 
and FR-I radio galaxies) is given by. 



In the case of Lacertids, relativistic beaming requires ^ ^ ^ l/F or equivalently 2r ^ 
^Lac ^ r. On the other hand, we suppose that off-axis counterparts verify 5 ^ 2/T (corre- 
sponding to an average angle value of 6* a; 60 deg for F > 1). Then, equation (22) rewrites 



This estimate can of course be complicated by an intrinsic luminosity distribution. It may 
be also that we cannot detect the unbeamed sources due to limited sensitivity of the instru- 
ment. However some other indicators such as the extended radio lobes power or the galaxy 
luminosity itself are not highly beamed, and can serve as an unbiased criterion to select 
samples. 

Capetti & Celotti (1999) studied a sample of 12 Lacertids and 5 FR-I sources with HST 
and compared the core luminosity ratio between objects sharing similar radiative properties. 
It clearly appears that the whole emission of the Lacertid cores is roughly 10^-10^ times 
brighter than the corresponding radio galaxy ones. Moreover Chiaberge et al. (2000) per- 
formed a similar work on a larger and more complete sample. They roughly obtained the 
same conclusions : the luminosity ratio between Lacertids and radio galaxy belongs to the 
interval lO^-^-lO^-^. Applying relation (23), we obtain typical values of F 2-5 for the bulk 
Lorentz factor. They also compare the broad band spectra of both classes of objects and 
they found that the spectra could be deduced by a simple Doppler boosting, but once again 
with modest values of the Doppler factor. 



The nearby giant eUiptical radio galaxy M87 (NGC4486, Zs ^ 0.00436) is the first (and 
for the time being unique) detected unbeamed radio-loud source at the TeV energy range. 




(22) 




(23) 



3.4. 



Detection of TeV unbeamed source — the case of M87 
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First detection was reported by the HEGRA collaboration with an integral flux above 250 
GeV at about 3.3% of the flux of the Crab Nebula (with a signiflcance of 4.7(7) during an 
high state (Aharonian et al. 2003; Bcilicke et al. 2004). Such TeV events are conflrmed by 
recent measurements of High Energy Stereoscopic System (HESS) (Beilicke et al. 2005). The 
powerful radio jet of M87 has been well studied in various wavelengthes from radio to X-rays, 
showing that the jet axis makes an angle between 30 and 40 deg with the line of sight. This 
angle is clearly large enough to ensure that the emission is unbeamed. Previous works based 
on the study of the proper motion of the VLBI knots (Biretta et al. 1995) or on the detailed 
analysis of HST and VLA observations (Lobanov et al. 2003) converge toward value for F of 
3-5 at the kiloparsec scale. The jet differential flux of a source expressed in the observer's 
frame can be written as function of the intrinsic differential luminosity as 

Mu;9,z)^{l + z)S' (24) 

with p = 1^*5/(1 + z) and where di{z) ~ zc/Hq is the usual luminosity distance. We now 
consider two different versions of the same intrinsic object, a beamed one corresponding to 
a blazar and an unbeamed one corresponding to a radio galaxy. In this case, the ratio TZ of 
the observed photon fluxes above some threshold frequency z/*'^'' writes. 



Lti.rdu[T^{iy)/hu] \zuj \1 + Zuj^ \Sby ^25) 



k{zb,Zu;a) 



where the index „ (rcsp. 5) refers to the observed unbeamed (rcsp. beamed) quantities, and 
where we suppose that the high energy spectrum can be expressed as a simple power-law 
with a photon index a. 

For beamed sources, the Doppler factor can be written as ^5 ^ 2F while for the unbeamed 
case one has 6u = 1/F(1 — /? cos 6') < F with /5 ~ 1 — 1/2F^. Finally we can express the bulk 
Lorentz factor as a function of 9 and the observational parameters only 

2(1- cos ^) > ■ ^^^^ 

A raw approximation of the previous expression is 



r(^) ^ I 



7^ ' 



showing the ~ 1/^ functional dependence of F and its slow power- law variation with TZ (or 
k). For instance, for a typical value of a = 2.5, a factor ten on TZ implies only a factor 
{10^^^ ^ 1.29) on F. 
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In 1997 April flaring period, the TeV blazar Mrk 501 {zg ^ 0.034) became roughly 8 
times as bright as the Crab Nebula as reported by the French collaboration CAT (Djannati- 
Atai' et al. 1999). Assuming M87 is an unbcamcd counterpart of Mrk 501 with an angle 
30deg ^ Om87 ^ 40deg we obtain 4 ^ F ^ 5.3. Again we find that the luminosity ratio is 
compatible with modest values of the Lorentz factor. Due to the increasing sensitivity of the 
present and the next generation of the Imaging Atmospheric Cerenkov Telescope Arrays, the 
detections of more and more TeV radio galaxies should help us to constrain the dynamics of 
the emitting plasma at the subparsec scale in a more reliable statistical way. 

3.5. Summciry 

All the above considerations show that observational data are compatible with the 
beaming model only if the bulk Lorentz factor for the X-ray and TeV emitting part of the 
object is relatively low, between 3 and 5. This value reproduces correctly the luminosity ratio 
and the statistical number of sources (which are a priori independant factors). Conversely, 
a value of F = 12.5 which is the minimum typical value derived from the one-zone modeling 
approach, would lead to a luminosity contrast oi uj ^ 10^-^-10^. This latter estimation is 
clearly not compatible with the previous observations, ascertaining the " Bulk Lorentz factor 
crisis of Te V blazars" . In the following, we will examine some suggestions made by various 
authors to solve the crisis. 



4. How to solve the crisis 

4.1. Two pattern model 

Chiabcrgc ct al. (2000) and Trussoni ct al. (2003) argue that a jet velocity structure 
can solve the problem of the BL Lac/FR-I unification scheme. They consider a (med- 
)relativistic external layer and a fast internal spine which dominates the emission in the case 
of a favorable alignment along the observer's line of sight, i.e. in the blazar case. Although 
similar in appearance to the two-flow model of Pelletier (1985) (see below for details), it 
differs by the fact that both flows are relativistic, one with a "low" Lorentz factor (around 
3) and one with a high Lorentz factor (at least 10). In the following, we consider the same 
approach considering a two- components modeling of the velocity structure, where a fast 
inner structure is supposed to be surrounded by a slow one. Each of these components is 
respectively characterized by a bulk Lorentz factor Ff and Fg. As we saw, the radiative 
emission of the moving source with a bulk Lorentz factor F is beamed in a cone sustained 
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by a solid angle SQ — tt/F^ along the motion. Therefore the emitted radiation appears to 
be Doppler boosted when the jet lies into Sfl around the observer's line of sight. In this 
case, the luminosity contrast between the two parents population (Lacertids and FR-I radio 
galaxies) writes, 

Tl^w^ (^,^f)^ (28) 

where right and left bound correspond to the case where fast velocity component respectively 
dominates or not the emission. Unification models are sensitive to the slow component only, 
so (Fs = 3 — 5) (Hardcastle et al. 2003; Trussoni et al. 2003; Chiabcrgc ct al. 2000; Urry & 
Padovani 1991). Therefore, assuming Fg ~ 4 we obtain 10^'^^ ^ w ^ 10^ if Ff = 25 and 
10^-«^ ^ ^ lO^-^ if Ff = 10. 

The results give a possible solution to the luminosity ratio problem, being more compatible 
with observations. But in the next section we will examine the consequences of such a 
velocity structure on the detection probability of TeV emitters among BL Lac objects. 



4.1.1. Statistics of detected sources 

Suppose a population E of sources randomly oriented per unit volume uq. Then the 
density number of sources oriented with an angle 9 — cos~^ /i according to the observer's line 
of sight is dn/d/i = no/2. We define BL Lac sources as those seen into the 6Qs — '^1^1 cone 
and therefore the TeV emitters as the part of Lacertids lying into the 50,^ — Tr/Ff . Therefore 
the probability of detecting a BL Lac object in E writes 

1 /"'^o dn 1 
VLs.^r{ii^lio) = - d/x— = — , (29) 
no Ji d/i AVj 

where Hq = 1 — l/2Fs. The probability that a BL Lac object a; e E is also a TeV emitter 
source {i.e. a source which the emission is dominated by the fast inner structure) is given 
by the conditional probability. 



T^TeV/Lac = V{uj G TeV|u; e Lac) = 



(30) 



The probability of detecting n TeV emitters among a population of N Lacertids is given by 
the usual binomial probability law, 

nn/N) = ^TeV/Lac(l " ^TeV/Lac)^-^ (31) 

and implicitly depends on the value of the ratio Fg/Ff. It is more convenient to express the 
probability of detection of at least n TeV emitters among the same sample of N Lacertids 
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which writes, 



k=N 



V{N ^no^n) = Y^ V{k/N). (32) 



k=n 



4- 1.2. Applications 

For a given value of {n,N), requiring that V{N ^ uq ^ n) is larger than an a pri- 
ori probabihty Vq constrains the space parameters (FajPf). The latter inequality leads to 
eliminate parameter's region lying above a straight line which corresponds equivalently to 
a constant value of Vq or of the ratio Fg/Ff (see cq. [32], [31] and [30]). Further restric- 
tion come from the 7— rays transparency argument and FR-I/Lacertids unification models 
as developed above. 

1. Firstly, the 7— rays transparency argument developed in the first part of this work 
directly constrains the value of the fast component as it requires a minimum value of 
the Doppler factor ^min, and therefore Ff ^ 5-^i^/2. We have shown that dy^in ~ 50 
exclude all part of the (Fg, Ff) lying above Ff = 25. 

2. Secondly, basic statistics argument based on the number FR-1 radio galaxies regarding 
the Lacertid one and the comparison of luminosity distribution of the previous popu- 
lations constrain the value of the slow component to reasonable values less than Fg 7 
(Urry & Padovani 1995; Chiaberge et al. 2000). 

We test this result on the catalog of BL Lac objects from Padovani & Giommi (1995). At 
Zs ^ 0.13 they report 29 Lacertids with known redshift. Setting {n — l,N — 29), Vq — 1% 
and recalling that Fg^max = 7 and Ff = 25, the intersection of all listed previous constraints 
reduces to null region (see figure 5). Even with the hypothesis of a structured flow, a large 
value of the Lorentz factor recquired by one-zone homogeneous models is clearly untenable 
(excluded with a confidence level of 99%). 

We demonstrate that even if the two-components velocity structure can give a satis- 
factory answer to the luminosity problem of the Lacertids even with large value of Doppler 
factor required by high energy emission models, it fails to explain the detection statistics of 
the TeV emitters among the BL Lac object population supposed to be off-axis FR-I sources. 
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5. Discussion 

5.1. Inhomogeneous models 

Altogether, the previous considerations lead to a serious paradox, where a high Lorentz 
factor larger than 20 seems mandatory to avoid strong 7 — 7 absorption, whereas all other 
facts tend to favour modest values around 3. The only way to solve the discrepancy seems 
to give up the implicit assumption of all one zone model, i.e. the fact that all photons are 
produced co-spatially and simultaneously in some characteristic region of size R. Alterna- 
tive to one-zone models have already been discussed in the literature. For instance, in the 
"blob-in-jet" model (Katarzyiiski et al. 2001, 2003), low energy photons are produced in a 
continuous jet and only the high energy ones are produced in a spherical blob. This allow to 
fit the overall spectrum with a smaller Doppler factor of around 15. Another possibility is to 
take explicitly the variability and use a time-dependant model to reproduce the data. Again, 
the constraints arising from 7 — 7 opacity can be somewhat released because soft photons 
are emitted at a later stage than high energy ones. As has been remarked by Ghisellini et 
al. (1985), time-dependant model will produce effects comparable to inhomogeneous ones. If 
the evolving source is moving at relativistic velocity, and that many flares are contributing 
to the emission, the overall system will be in fact an stratified jet composed with many 
"one- zone" regions in a different evolutionary stage. However, none of these models do use 
bulk Lorentz factors as low as 3. 

We are led thus to consider models where photons are distributed along a jet in a 
continuous structure, instead of filling a spherical source. In this case, the luminosity is 
proportional to the photon density times the lateral surface of the jet, which is 2TrRjhj = 
27r AR'j, where Rj and hj are the typical jet radius and length at the emission region, and 
A — hj/Rj is an aspect ratio of the source. For a self similar jet for which all quantities 
(radius, magnetic field, particle density etc..) are described by power law as a function of 
the distance z, one expects hj ~ Zj, where Zj is the distance of the emitting region from the 
center. It follows that A ~ Zj/Rj ~ 9j^, where 9j is the typical opening angle of the jet. One 
can see that for a given synchrotron luminosity and photon density (implying the same IC 
luminosity), one must conserve the quantity AR^, so the typical radius of the jet, and hence 
the 7 — 7 optical depth will be reduced by a factor A^^"^ with respect to a spherical source. 
This simple geometrical modification helps thus to increase luminosity without increasing 
optical depth. Furthermore, the particle distribution needs not to be the same all along the 
jet. Rather one expects a gradual cooling of the particles, the overall spectrum being the 
envelope of all slices of the jet. The local photon spectrum can thus be different from the 
observed one, and particularly the local soft photon density can be much lower, helping again 
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to reduce the 7 — 7 optical depth. As we shall see, all these factors can offer a clue to the 
Bulk Lorentz factor crisis, but imply strong constraints on the physical picture of relativistic 
jets. 

5.2. Theoretical implications 

5.2.1. Local photon density 

Considering the above constraints, we will take the opposite attitude, considering that 
the value of the bulk Lorentz factor is constrained by the unification models and the detec- 
tion of unbeamed sources to be around 3. The typical high energy emission zone as defined 
above is equivalent to the superposition of A spherical sources with individual luminosities 
L,^/A ~ OjL,^. Therefore, all previous equations in section 2 are still valid provided we re- 
place the observed luminosity L^, by 9jLi,. Using equation (12), we conclude that all opacity 

—1/2 

constraints remain unchanged if we replace the optical depth Tjj by 9^ Tjj, which is of 
course in accordance with the estimate made in the previous paragraph. So we can use the 
figure 3 with slightly different values of r^^. The typical angle 6j must be of the order of 
10~^ — 10~^, so the optical depth will be reduced by a factor between 3 and 10. In the 

— 1/2 

following, we will still use the same line 6j r^-y a; 1 to constrain the optical depth, meaning 
that r^^ ^ 0.1 to 0.3. 



5.2.2. Local photon spectrum 

We can thus put an upper limit on the soft photon luminosity corresponding to this 
value and a Doppler factor of 3, which constrains the soft photon luminosity at an energy 
Et = S'^/Ec ~ 10 eV. As the spectrum is by definition approximately the same in all the char- 
acteristic emission region, we can thus estimate the local photon spectrum by interpolating 
between the peak synchrotron luminosity and the above upper limit. Inspection of figure 
3 shows that the spectral index between 10 eV and 100 keV is very close to 1/3, which is 
characteristic of a quasi monoenergetic distribution; an example of such distribution is pro- 
vided by the quasi-maxwellian or "pileup" distribution (Henri & Pelletier 1991; Schlickeiser 
1985; Sauge & Henri 2004), which is a natural outcome of some acceleration processes like 
second order Fermi acceleration or magnetic reconnexion. This distribution is not the usual 
power law often claimed to exist in AGN, and which is naturally produced in MHD shocks. 
Rather than localized shocks, the assumption of low Lorentz factor leads to a picture of a 
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continuous jet filled by relativistic particles, continuously reheated by a diffuse acceleration 
mechanism. 



5.2.3. Pair production 
The local synchrotron spectrum can not be harder than the monoenergetic one; so figure 

I/O 

3 proves that this implies a lower limit to the quantity 9^ r-y-y > 1. Thus the limit on 7 — 7 
optical depth can not be very low, unless we have an extremely well coUimated jet that 
is not supported by the general FR-I morphology. Modest collimation factors imply that 
'^77 ^ O-l- This supports the formation of an electron-positron pair plasma in the acceleration 
site. If the acceleration is not localized, which is suggested by the picture of a continuous jet 
filled by a pile-up distribution, the pairs created by 7 — 7 interaction can not avoid being 
reaccelerated and will trigger a pair cascade (Henri & Pelletier 1991). These constraints are 
thus suggestive of a inner continuous pair dominated, jet-like emission zone, maintained at 
a relativistic temperature. 

5.3. Compatibility with the two-flow model 

All the previous considerations find a natural explanation in the context of the two-fiow 
model, which was proposed to account for the formation of relativistic jets in AGN : in this 
model, extragalactic jets are in fact the results of a double structure: a first jet, not highly 
but only mildly relativistic {v 0.5c), is emitted by a MHD mechanism by a large scale 
magnetic field anchored in an accretion disk (Blandford & Payne 1982; Ferreira & Pelletier 
1993a,b, 1995) ; this powerful, but weekly dissipative jet, can sustain a MHD turbulence able 
to accelerate non thermal particles. These particles will produce synchrotron and gamma- 
rays photons, and if the optical depth becomes large enough, these photons will trigger an 
intense pair cascade leading to a dense pair plasma in the empty "throat" of the jet. We 
have shown in previous works that this pair plasma will be spontaneously accelerated to 
relativistic velocities even if the surrounding jet is not highly relativistic by itself, by the so 
called " Compton Rocket effect" , which is a recoil effect associated with anisotropic IC pro- 
cess originally introduced by Odell (1981). The Compton Rocket effect has been shown to 
be inefficient to accelerate an isolated relativistic plasma because the cooling time is always 
shorter than the bulk acceleration time (Phinney 1982). In the two-fiow model however, the 
heating by the surrounding jets compensates for the cooling and the pair plasma remains 
relativistic over large distances (Marcowith, Henri, & Pelletier 1995). 
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Detailed calculations of the Compton Rocket effect in this configuration (Renaud & 
Henri 1998) show that the pair plasma accelerates gradually in the vicinity of an accretion 
disk, being maintained to a quasi equilibrium Lorentz factor Fbeq ~ (z/rj)^/^, where is 
the inner radius of the accretion disk (3 gravitational radii for a Schwarschild black hole). 
The equilibrium Lorentz factor is defined by the fact that the photon field of the accretion 
disk, seen in the comoving frame, appears to be nearly isotropic due to relativistic aberra- 
tion. It grows slowly with the distance,because the field becomes more and more anisotropic. 
The acceleration continues until the photon density becomes to low to efficiently accelerate 
the plasma. Then the plasma decouples from the ambient radiation field and ends up with 
an asymptotic balistic motion at constant Fb — > Fboo, which depends on the disk lumi- 
nosity and the particle energy distribution. For a relativistic energy distribution function 
n{j) oc 7^exp(— 7/7), the asymptotic bulk Lorentz factor is approximately Fboo ~ i^-slY^'^ , 
where ig = LsaT^ii/ATrrrigC^ri is the soft photon disk compactness and 7 is the characteristic 
energy of the pileup depending on the details of the acceleration/cooling processes (Renaud 
& Henri 1998). 

The first interesting feature in this model is that it predicts naturally a gradual accel- 
eration from the core. The value of F^ ~ 3 is naturally obtained at ~ 100 Vg, which is a 
typical distance where gamma-ray emission seems to occur, based on variability arguments. 
Thus low Lorentz factor are not surprising in this model, but are explained naturally. As a 
matter of fact, very high values of 20 near the core would be difficult to explain in this frame! 

The second one is that the asymptotic bulk Lorentz factor is controlled by the density 
of the photon field emitted by the accretion disk. For BL Lac objects and FR-I galaxies, the 
disk luminosity is known to be much lower than luminous FSRQ and FR-H galaxies, by a 
factor around 10~^. One would expect thus a lower asymptotic Lorentz factor for BL Lac 
object in average, which would help to understand the absence of superluminal motion in 
TeV blazars. As a matter of fact, numerical estimates show that the expected asymptotic 
Lorentz factors are between 10 and 20 for near-Eddington accreting supermassive black holes, 
whereas they are rather between 5 and 10 for low luminosity AGN. We note that bulk Lorentz 
factors around 5 are indeed observed in M87, which would mean that the decoupling occurs 
at some thousands Schwarzschild radii from the core. Unification models are compatible 
with slowly accelerating jets, the inner (X-ray emitting) jets having bulk Lorentz factors 
around 3 and the outer radio jet having a larger Lorentz factor around 7 (Urry & Padovani 
(1995)). Again this is perfectly compatible with the predictions of the two-fiow model, with 
an inner jet emitting X-ray and TeV radiation with a modest bulk Lorentz factor, and an 
outer jet responsible for radio emission with a higher one. Also we note that there is no 
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need for deceleration to explain FR I mildly relativistic jets : even if the large scale jet has 
only a moderately relativistic velocity v ~ 0.5c, this can be attributed to the "slow" MHD 
component surrounding the relativistic beam, the latter being dissipated at kpc scale. 

An inhomogeneous model offers also a convenient explanation for the lack of obvious 
correlation between X-rays and gamma-rays variability. If the magnetic field is varying 
along the jet, the photons with a given energy could be produced by electrons with different 
energies and locations. If several flares contributes to the observed spectrum — which is 
necessary to account for the global spectral shape in case of a monoenergetic distribution 
— a complicated variability pattern could emerge. This is much less easy to understand in 
the homogeneous steady state models. Thus we think that inhomogeneous models, although 
more complicated to compute, seem to be unavoidable to explain the spectral and temporal 
features of TeV blazars' emission. 



6. Conclusion 

We have investigated in detail the so-called " Bulk Lorentz factor" crisis of TeV blazars, 
which seem to imply an incompatibility between a high Lorentz factor required to insure 
gamma-ray transparency, and a low Lorentz factor deduced from statistical arguments and 
luminosity contrast, including the detection of the non-blazar TeV source M87. We show 
that the transparency argument is common to all one-zone models, and that the only way 
of solving the paradox is to consider inhomogeneous jet models, where all photons are not 
produced cospatially. The spectrum is then the spatial convolution of different jet slices, 
and the opacity problem can be avoided by invoking geometrical arguments and harder 
local photon spectrum. We show however that for modest values of geometrical beaming of 
the jet, which seem natural considering the morphology of FR-I galaxies, the optical depth 
for 7 — 7 absorption can not be very low, even for a local quasi-monoenergetic particle 
distribution. This has profound implications on the physics of the jet : the acceleration 
mechanism must be distributed all along the jet, and is more probably insured by second 
order Fermi mechanism or reconnexion sites than by localized shocks . A moderately high 
value of 7 — 7 optical depth implies a fair production rate of electron-positron pairs, which 
are likely to be reacceleratcd by the acceleration process to trigger a pair cascade. All this 
features are natural consequences of the two-flow model, which attributes the relativistic 
phenomena (high energy emission and supcrluminal motion) to the formation of such a pair 
plasma inside a powerful, but mildly relativistic jet insuring the confinement and the heating 
of the relativistic beam. The bulk Lorentz factor is also well in accordance with a continuous 
acceleration along the jet by the Compton Rocket effect, which predicts naturally Fb ~ 3 
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at a hundred Schwarzschild radii from the core. We conclude that all observational facts 
are more in accordance with light, moderately relativistic leptonic beams than with highly 
relativistic baryonic jets. 

Remarks of an anonymous referee helped to improve significantly the final version of 
this paper. L.S. would like to thank all members of the IPNL team of the SNFactory 
collaboration. Part of the simulations reported here has been performed at the "Centre de 
Calcul Intensif de I'Observatoire de Grenoble". 
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Table 1. Mkn 501 1997 April- 16 observationnals parameters 



DiRB status 


Synchrotron bump 


Inverse Compton bump 






logio/^r^ logiof^.^) 

\ " / max 


\ ^ / max 


'"max 


reddened 
de-reddened 


19.3 ±0.1 -9.11 ±0.04 


26.4 ±0.1 -9.46 ±0.04 
26.6 ±0.1 -9.04 ±0.04 


0.4467 
1.175 
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Table 2. Values of function a^^{P) and the modified function modified function a-y^{P) as 

function as the spectral index (5 in uFi,. 



p 


1 
2 


1 


4 
3 


a^^{(3) 0.122 


0.236 


0.583 


1.397 




0.043 


0.583 


1.613 
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photons "created" by pairs with 
the same Lorentz factor 




Photons involved in the same 
pair creation procees 

Fig. 1. — Working figure - High energy photons with reduced energy e* interact preferentialy 
with soft photons of energy 1/e* to create new pairs. In the Klein-Nishina scatering regime, 
an ultrarelativistic particle of reduced energy 7 = 5* "create" a high energy photon with 
the same energy e*, and in the same time, other pairs with the same energy can create soft 
photons by synchroton process of energy e* = {3eh/47rmlc^) Bj"^ = [B/Bq) 7^. 
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Fig. 2. — Spectral energy distribution of Mrk 501 during the flaring period in 1997 April, 
showing the simultaneous data taken by the Beppo-SAX instrument (Pian et al. 1998) and 
by the CAT imaging Atmospheric Cerenkov Telescope (Djannati-Atai et al. 1999; Barrau et 
al. 1998). About high energy data points, filled gray circles are the CAT observed ones while 
open squares are unabsorbed ones, corrected from our estimation of the DIrB attenuation 
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Fig. 3. — Constraints on the local shape of the synchrotron spectrum of Mrk 501 during 
the high 1997 April 16 high state. Grayed polygon is obtained considering the gamma 
transparency argument. It is defined by the zone where opacity lying into the interval 
T-y-y X y4^/^ G [0.1, 1] where A is the aspect ratio. For an homogeneous spherical blob A = 1 
while in the case of a jet A = l/Oj where 6j is the characteristic opening angle of the jet. 
Constraint coming from the typical variability time scale leads to the most left straight thick 
line. Also represented in dotted-dashed line, a spectrum with a spectral index equals to 
4/3 in vFjj resulting from the emission of a (quasi-)monoenergetic distribution of electrons 
and/or positrons. Left (resp. right) pannel correspond to the de-reddened (resp. reddened) 
case. 
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FR-I [is 




(spine) 



Fig. 4. — Sketch of the Two pattern modeL It considers a relativistic external layer charac- 
terized by a Lorentz factor Fg and a fast internal spine with Fg > Ff . Considering relativistic 
Doppler beaming, an object viewed with 9 > l/Fg refers to a FR-I radio galaxy. Conversely, 
if ^ l/Fg the source is seen as a BL Lac object and more precisely, if ^ 1/Ff the fast 
inner component dominate the emission. In this latter case one deals with a TeV blazar. 
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Fig. 5. — Left panel. Space parameters (Fg, Ff) constrained by the statistics of the number 
n of TeV blazars in a given population of N observed Lacertids. Each Une corresponds 
equivelently to a constant value of Vq = V{29 ^ rio ^ 7) or of the ratio Fg/Ff. Here 
represented Vo — 95%, 90%, 50%, 1% and 0.1% (see text for more details). Right panel. Same 
as previous panel but combined with the constraints coming from 7— rays transparency 
argument and FR-I/Lacertids unification models. The first one ehminates all the region 
lying above the Ff min = 25 while the seconde one suppress the right part of the parameters 
space Fs ^ Fg^max = 7. In this case, the allowed region compatible with a probability of 
detection of at least 7 TeV blazars in a population 29 Lacertids is strongly improbable. 



